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Using A Multi-Phase Transport (AMPT) model, we study J/tp production from interactions 
between charm and anti-charm quarks in initial parton phase and between D and D mesons in 
final hadron phase of relativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC). 
Including also the inverse reactions of J/ip absorption by gluons and light mesons, we find that 
the net number of J/tp from the parton and hadron phases is smaller than that expected from the 
superposition of initial nucleon-nucleon collisions, contrary to the J/tp enhancement predicted by 
the kinetic formation model. The production of J/ip is further suppressed if one includes the color 
screening effect in the parton phase. We have also studied the dependence of J/ip production on 
^— J ■ the charm quark mass and the effective charm meson mass. 
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A central interest in modern nuclear physics is to produce and study the Quark-Gluon Plasma (QGP) predicted 
\ by the Quantum Chromodynamics (QCD) [Q. Studying the properties of the QGP and the deconfinement phase 
• transition not only is important for understanding the QCD but also has astronomical implications Q. Experiments 
' at the Super Proton Synchrotron (SPS) involving collisions between heavy nuclei have shown that a large amount of 
, energy is deposited around mid-rapidity, and the energy density may already be sufRciently large for the formation 
T— I of the QGP during the initial stage of collisions More recent experiments at RHIC at Brookhaven National 
Laboratory, where the collision energy is much higher than at SPS, have also shown possible effects due to the 
formation of a partonic matter Q]. 

To find the signal for the Quark-Gluon Plasma in relativistic heavy ion collisions, Matsui and Satz proposed to 
study J/4> production in these collisions Based on the fact that the effective potential between charm and 
anti-charm quarks changes in the QGP due to color screening effect, no bound states can be formed between them 
beyond a certain critical temperature, which is somewhat higher than the deconfinement temperature. As a result, 
^ ■ J/ip production is expected to be suppressed if the QGP is formed in heavy ion collisions. There are other possible 
, mechanisms for J/ip suppression 1^-^ in heavy ion collisions, e.g., J/?/; can be destroyed by collisions with incoming 
nucleons or with gluons in the initial partonic matter |P,|lo| or with comoving hadrons in the hadronic matter . 
The observed abnormal suppression of JM in central Pb-|-Pb collisions at SPS may require the formation of the 
■ Quark-Gluon Plasma in these collisions [|l§,|| . 
5-H ' For heavy ion collisions at RHIC energies, unlike in previous fixed target experiments at SPS, multiple pairs of 
charm-anti-charm quarks can be produced in one event, and these charm and anti-charm quarks can interact and 
produce the J/tp . This mechanism is the inverse reaction of J/ ip absorption by gluons. Using a kinetic approach, 
it was found that despite J/ip dissociation due to color screening, this reaction is important in the quark-gluon 
plasma with temperatures between those for J/tp dissociation and the deconfinement, and would lead to an enhanced 
production of J/ip in heavy ion collisions at RHIC. A possible signal for QGP at RHIC has thus been changed from 



|23[ , which has been extensively used 
icts a comparable J/ip production. The 



J/4' suppression to J/tjj enhancement. The statistical fragmentation model 
for studying particle production in heavy ion collisions at SPS energies, prec 
main assumption of this model is that J/ip formed during hadronization of the QGP is in chemical equilibrium with 
charm mesons. 

In the kinetic formation model, the heavy ion collision dynamics is treated schematically. For a more quantitative 
study of this new J/ip production mechanism, we study J/ip production from both the charm-anti-charm quark 
interactions and the charm-anti-charm meson interactions using the AMPT model [p4|-p6|. In particular, we consider 
central (b=0) Au+Au collisions at the highest RHIC energy ^/s — 200 AGeV. With this dynamic transport model, 
we find that the net number of produced J/ip from the parton and hadron phases is smaller than that expected from 
initial nucleon-nucleon collisions, contrary to the J/ip enhancement predicted by the kinetic formation model. The 
production of J/ip is further suppressed if one includes the color screening effect in the parton phase. 

In Section ||, we give a short description of the multi-phase transport model (AMPT). Results are presented in 
Section |l| for J/ip production from the initial parton phase and the final hadron phase. In this section, we also 
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discuss the dependence of J /ip production in the parton phase on the charm quark mass and in the hadron phase on 
the charm meson mass. Finahy, we summarize the paper in Section IV. 



II. MULTI-PHASE TRANSPORT MODEL AND J/^j PRODUCTION MECHANISMS 

The muhi-phase transport model p^p6[ that we use for the present study is a hybrid model based on three Monte- 
Carlo models for the three stages of relativistic heavy ion collisions, i.e., the Heavy Ion Jet Interaction Generator 
(HIJING) p7| , p8[ for the initial conditions, Zhang's Parton Cascade (ZPC) for the parton evolution, and A Rel- 
ativistic Transport (ART) model |]30| , ^ for the hadron evolution. In the HIJING model, nucleons are distributed 
according to the Woods-Saxon distribution. The parton distribution in a nucleon is adjusted according to the trans- 
verse position of the nucleon and a shadowing function based on the Mueller-Qiu parton recombination mechanism 
p2| . These partons undergo hard or semi- hard collisions and generate an initial parton system after the passage 
of two colliding nuclei. HIJING provides the momentum space information of the partons produced in the hard or 
semi-hard collisions. The formation time of a parton is generated according to the Gyulassy-Wang model of parton 
coherent production . The position of a physical parton is then obtained from the position of its parent nucleon by 
adding a displacement given by its velocity times the formation time. The resulting parton phase space distribution 
is used as the initial conditions for starting the ZPC parton evolution. In the present version of the ZPC model, 
only elastic gluon-giuon scattering arc included with a cross section of 3 mb, which is consistent with the estimated 
average screening mass in the partonic matter formed in the collisions. When the partonic matter stops interacting, 
partons are connected back with their parent nucleons, which also suffer soft interactions through string excitations. 
The Lund fragmentation model is then called from the HIJING model to convert these wounded nucleons into 
hadrons. These hadrons become physical after a formation time given by the average freeze-out proper time of partons 
in a string plus an additional 0.7 fm/c to account for the time that is required for producing quark-anti-quark pairs. 
The ensuing hadron evolution is modeled by the ART model that takes into account interactions among nucleons, 
anti-nucleons, mesons, and resonances. For a better description of experimental hadron spectra, the original ART 
model has been improved by including the K* meson interactions, the anti-nucleon production, and the p meson 
interactions. 

We have previously used this multi-phase transport model to study the suppression of primarily produced J /ip in 
heavy ion collisions at RHIC These J/tp^s are produced in nucleon- nucleon collisions during the initial passage of 
two colliding nuclei. Because of color screening, these J/ip^s are not physical ones but correlated pairs of charm-anti- 
charm quarks close in phase space. They become physical J/tJj^s if they remain close when the partonic matter cools 
below the J/^p dissociation temperature. Otherwise, they will combine with nearby light quarks to form charm mesons 
at hadronization. Gluons can also destroy the phase space correlation between the charm quark and anti-quark pairs 
and hence destroy the J/ip. It was found that finite space size and lifetime have significant effects on J/ip suppression, 
and different mechanisms lead to different suppression factors. Also, the dependence of J/ip suppression on the size of 
colliding nuclei was shown to be useful for differentiating the different suppression mechanisms. Similar results have 
also been obtained in other studies ||3^,0 . 

In the present study, we include the reactions of J/^p production from charm- ant i- charm quark interactions in the 
parton phase and from charm meson interactions in the hadron phase. As the primary J/^p number is much smaller 
than the number of other particles participating in the equilibration, primary J/ip evolution won't affect secondary 
J/ip production. We do not include primary J/tp^s in the following calculations. 

The charm and anti-charm quarks are generated from the PYTHIA model , with a cross section of 350 /.ih for 
producing a pair of charm-anti-charm pair in a nucleon-nucleon collision. This cross section does not include nuclear 
shadowing effect of gluons, so the present study gives an upper bound for charm meson and J/tp production in heavy 
ion collisions. In the partonic stage of heavy ion collisions, J/ip's are then produced from charm-anti-charm quark 
interactions through the reaction c -I- c — > J/ip + g. The J/ip can also be destroyed by gluons through the reaction 
g + J/ip — !■ c -|- c. We further include the dissociation of J/ip due to color screening in the partonic matter. To treat 
this process, we follow the method used in Ref. ]3^ by introducing a critical radius in the transverse plane of a heavy 
ion collision, within which the parton energy density is larger than the critical energy density for J/ip dissociation 
by color screening and the J /i]^ thus cannot be formed. The critical radius decreases as the partonic matter expands. 
The time evolution of this critical radius in Au+ Au collisions at -^s = 200 AGeV has been parameterized in Ref. . 
In addition to producing the J /ip, charm and anti-charm quarks also undergo elastic scattering with gluons. This 
elastic cross section is taken to be 3 mb in the present study. 

After partons stop interacting, charm and anti-charm quarks are converted to D and D mesons. This conversion 
is carried out using a delta function fragmentation scheme. In this approach, the momentum of a charm meson is 
taken to be the same as the momentum of its parent charm quark. A formation time of 1 fm/c is also introduced in 
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combining the charm quark with a nearby Ught quark. Because of the large mass of charm quark, the delta function 
fragmentation scheme gives a reasonable description of charm quark fragmentation in hadron reactions. For simplicity, 
we do not differentiate between D mesons and D* mesons. Instead, we vary the D meson mass to study the influence 
of different production thresholds for J /ijj production from charmed hadrons. 

In the hadronic phase, D and D mesons can interact to produce a J /if) and a light meson through the reaction 
DD J /ipM , where M denotes a light meson. The J /ip can also be destroyed by light mesons through the inverse 
reaction J/tpM —^ DD. We only include reactions that involve tt mesons or p mesons as they are the most abundant 
ones in the nearly baryon free matter formed in heavy ion collisions at RHIC. The cross section for light meson 
destruction of J/^ is taken to be 3 mb as in previous transport study at SPS energies [p"3[-p^. This cross section 
is consistent with those from recent theoretical studies based on effective hadronic Lagrangians [^9|-^, the quark- 
exchange model and the QCD sum rules The cross section for the inverse reaction of J/t/j production from 
D mesons is calculated according to the detailed balance relation. As we do not differentiate between D mesons and 
D* mesons, the spin and isospin degrees of freedom for both D and D* mesons are used in evaluating the inverse 
reaction cross section. This is equivalent to including all possible combination of D and D* states. We will vary the 
D meson mass to study the effect due to changing J/ip production threshold. 

In earlier studies based on kinetic models p4|, JM production from the hadron phase of relativistic heavy ion 
collisions was found to be important only at LHC but not at RHIC as the cross section for DD — > J/ipM used 

in these studies was smaller than that used in present study. 

Although multiple pairs of charm quarks are expected to be produced in each Au+Au collision, the number of 
charm quarks is small compared with other particles. For example, the rapidity density of charm quark pairs is 
about 1.5 at = 200 AGeV. The expected number of J/tj} per event is even much smaller. To obtain sufficient 
statistics for J production, charm particles and the J /ip are treated in the transport model by the perturbative 
method @,|^, i.e., neglecting their effects on heavy ion collision dynamics. This method has been used extensively 
for studying rare particle production in heavy ion collisions at low energies. Explicitly, we include many charm quark 
events from primary nucleon-nucleon interactions in the evolution of partons from a single HIJING event. However, 
only charm quarks from the same event are allowed to interact among themselves, and the usual method of treating 
binary collisions in transport models is used in treating these collisions. For collisions between charm and uncharm 
particles, we keep their effects on the charm particles but neglect those on the uncharm particles. In this way, the 
dominant computation time for the interactions between uncharm particles remains the same as for a single charm 
event. Since there are much more uncharm particles than charm particles, e.g., about 300 gluons verse about 3 charm 
and anti-charm quarks per unit rapidity, neglecting the effects of charm particles on the heavy ion collision dynamics 
is expected to be a very good approximation. 



III. RESULTS 



A. J/tp production in parton phase 

We first study J /ipi production in the parton phase. Unless otherwise indicated, the rates and numbers in the 
following are all rapidity densities averaged with \y\ < 1. Fig. |l| gives the time evolution of the production and 
destruction rates per unit rapidity for J/ip's with < 1. Also shown is the rate per unit rapidity for collisions 

between charm and anti-charm quarks which would have produced the J/ip in the absence of color screening effect. 
These unsuccessful collisions mainly occur before 0.5 fm/c when most partonic matter have effective temperatures 
above the J/ip dissociation temperature. The screening effect largely diminishes after 1 fm/c when the parton density 
decreases due to expansion. At this time, the J/iJj production rate starts to increase quickly, and this is followed by 
an increase of the J/ip destruction rate. After reaching their peak values, both the production and destruction rates 
decrease and become comparable at about 4 fm/c. 

To see more clearly the relative importance of J ftp production and destruction in the parton phase, we show in Fig. ^ 
the time evolution of the number of produced J/ip and the number of destructed J/tp as well as the net J/ip number. 
The number of J/ 4' saturates at about 3 fm/c. It decrease slowly afterwards as the production and destruction rates 
become comparable as shown in Fig. |^. 

In the absence of color screening effect, collisions between charm and anti-charm quarks in the partonic matter that 
is above the J/tp dissociation temperature can also produce the J/ip. This leads to a very high J/ip production rate 
before 0.5 fm/c as shown in Fig. ^. The J/ip destruction rate during this time also reaches its peak value of about 
1/3 of the maximum production rate. In contrast to the case including color screening effect, the production and 
destruction rates become comparable already at about 0.5 fm/c. As shown in Fig. ^, these result in a much faster 
increase in the numbers of produced and destructed J/ip, and an earlier peaking in the net J/ip number than in the 
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case with color screening effect. Since the destruction rate after 0.5 fm/c is still slightly larger than the production 
rate, the net J/V' number slowly decreases from its maximum value of about 0.003 to a final value of about 0.0014 per 
event. This number is about a factor of two larger than the final net J /ip number for the case with color screening 
effect, which is about 0.0007 per event. 

The number of J /ip expected from primary nuclcon-nucleon collisions in Au+Au central collisions can be estimated 
from the production cross section in pp collisions. The cross section per unit rapidity at ^/s = 200 AGeV is 
0.63 /lb Eq]. Using the nuclear overlap function TAu+Au{b = 0) = 30 mb~^, we expect 0.019 J/il) per unit rapidity 
to be produced in Au+Au central collisions. This is much larger than the predictions of 0.0007 and 0.0014 from the 
dynamical multi-phase transport model. Our study thus predicts a suppression of J / ip production in the parton phase 
of heavy ion collisions at RHIC, and this is opposite to the enhancement predicted by the kinetic formation model. 



B. J/tp production in hadron phase 

In hadronic matter, the threshold for J/tp production is reduced in comparison with that in the partonic matter as 
the D meson mass is larger than the charm quark mass. This threshold effect is expected to compensate for the lower 
charm density in the hadron phase than in the parton phase. To demonstrate this threshold effect, we vary the D 
meson mass in studying J/ip production in the hadronic matter. Fig. |^ shows the time evolution of J/iJj production 
and destruction rates per unit rapidity in the hadron phase for different D meson masses. These results are obtained 
with the color screening effect included in the initial parton phase. The upper panel is for mo — 1-70 GeV, which is 
smaller than the D meson mass in vacuum. Wc use this smaller mass because theoretical studies have shown that the 
D meson mass may be reduced in medium |49 -52j. The middle panel is for mo = 1.87 GeV, which is the D meson 



mass in vacuum. Here, we treat D* meson mass as D meson mass, and the results are thus a lower bound for J/^ 
production. The lower panel is for mo = 2.01 GeV, which is the D* mass in vacuum. Treating D meson mass as D* 
meson mass in this case then gives an upper bound for J/ ip production. The destruction rates for all three D meson 
masses are seen to have similar magnitude. For the production rate, it is larger for larger D meson masses. 

In Fig. ^, we show the time evolution of the numbers of produced and destructed J/ip together with that of the 
net J/ip number. For mo = 1.87 GeV and mo = 2.01 GeV, the number of produced J/ip is always greater than 
the number of destructed J ftp, leading to a net production of J/^ from the hadron phase of relativistic heavy ion 
collisions. With increasing D meson mass, the destructed J/^p number increases slightly while the produced J/^p 
increases drastically. For mo = 1.70 GeV, the produced and destructed J/^ numbers are about equal, and the final 
net J/ip number is thus about the same as that from the parton phase. On the other hand, for mo = 2.01 GeV, the 
final J/tp number equals approximately the number of produced J/^p from the hadron phase. The final J/V' number 
reflects essentially the one produced from the hadron phase. 

Quantitatively, the final net number of J/ip per unit rapidity averaged over | < 1 is about 0.0007 per event for 
mo = 1.70 GeV and increases to 0.0019 for mo — 1-87 GeV and to 0.0053 for mo = 2.01 GeV. The time evolution 
of the net J/ip number in the three cases are similar as it is largely determined by the heavy ion collision dynamics. 
The J/tp number saturates at about 15-20 fm/c, which is about the lifetime of the hadronic matter. 

The results for neglecting the color screening effect in the initial parton phase are shown in Fig. [t] for time evolution 
of the J/V' production and destruction rates and in Fig. ^ for the time evolution of the numbers of produced and 
destructed J/^p as well as the net J/ip number. Although there are more J/ip^s destructed in the hadronic matter, 
particularly for mo = 1-70 GeV where the number of destructed J/^ is even larger than the number of produced 
J/ip in the hadron phase, the final J/ip number is larger than in the case with color screening in the parton phase. 
Their numbers are 0.0010, 0.0024, and 0.0057 for mo = 1.70, 1.87, and 2.01 GeV, respectively 



C. discussions 

The relative importance of J/ip production from the parton and hadron phases is more clearly seen in Fig. ^, where 
we show the number of produced J/ip per unit rapidity with |yj/^| < 1 as a function of the D meson mass. The 
results with color screening in the parton phase are shown by open circles for production from the parton phase and 
filled circles for the final J/tp number including production from the hadron phase. Without color screening in the 
parton phase, the corresponding results are given by diamonds. As expected, including color screening effect leads to 
a decrease in the final J/^p number. As the D meson mass increases, the final J/ip yield increases strongly due to 
production from the hadron phase, and the effect due to color screening becomes less important. The results using 
mo = 2.01 GeV give the upper bound for J/tp production while those using mo = 1-87 GeV give a lower bound. We 
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note again that the net effect of a reduced D meson mass of mn = 1.70 GeV is the destruction of J/'0's produced 
from the parton phase in the absence of color screening. 

Part of increased production of J /ip in midrapidity when the effective D meson mass increases is due to increase 
in the charm meson rapidity density during hadronization. This is shown in Fig. where we give the midrapidity 
densities of charm quarks before (lower dashed line) and after (lower solid line) parton evolution as well as of charm 
mesons before (upper dashed line) and after (upper solid line) hadron evolution. During parton evolution, the rapidity 
density of charm and anti-charm quarks decreases slightly due to collisions with gluons. Hadronization increases the 
midrapidity density of resulting charm mesons, and its effect increases with increasing charm meson effective mass. 
The midrapidity density of D and D mesons is seen to increase by 8% when mo increases from 1.70 GeV to 2.01 
GeV. However, this would lead to at most 20% increase in the JN) midrapidity density and cannot account for the 
more than factor of 5 increase in J/ "0 production shown in Fig. S. The dominant effect for such a larger increase is 
thus due to the reduced threshold for J /ijj production when D meson mass increases. 

The charm quark mass strongly influences the number of J/t/; produced from the parton phase. Fig. |l^ shows the 
J/V' rapidity density as a function of the charm quark mass. To get the upper limit of J/ij) production, the D meson 
mass is set to be 2.01 GeV and color screening is not included. The J /ip number increases by a factor of 72% when 
the charm quark mass increases from the value 1.35 GeV used in PYTHIA to 1.5 GeV used in the constituent quark 
model. If the charm quark mass is taken to be 1.8 GeV, we observe about a factor of 4 increase in J/ij) production 
from the parton phase compared with the rric = 1.35 GeV case. 

In the kinetic model approach of Ref. |2^, the midrapidity density of charm quarks in Ref. is about 2.5 because 
the largest rapidity range for about 10 pairs of charm and anti-charm quarks is 4 in that study. This is larger than 
that shown in Fig. |l^, which is about 3.48/2=1.74. Since J/ip production is roughly proportional to the square of 
the number of charm quarks, this can lead to about a factor of 2 difference between our results and those from Ref. 

. On the other hand, the cross section we use for J/ip production from the charm and anti-charm quarks is about 
a factor two larger than the peak value used in the kinetic model. The effects due to different midrapidity charm 
quark densities and J/V" production cross sections in these two studies are thus largely cancelled out. The different 
conclusions of J /tp suppression in our study and J /ip enhancement in the kinetic approach are results of different 
heavy ion collision dynamics used in the two studies. In our approach, the parton phase is generated from the mini-jet 
gluons from initial hard and semi-hard collisions between nucleons and is evolved using the transport model. In the 
kinetic approach, an equilibrated quark-gluon plasma is assumed to be formed in the initial stage and later expands 
according to the 1-D Bjorken hydrodynamical model p3| . Further studies are required to determine the validity of 
these predictions. 

IV. SUMMARY 

In this paper, we have studied the effects due to interactions between charm and anti-charm quarks in the initial 
parton phase and between charm mesons in the final hadron phase on J /ip production in central Au-|-Au collisions at 
RHIC. Using a multi-phase transport model, we find that because of the high density of charm quarks at early times, 
the J /ip can be produced from the interactions of charm quarks. This is consistent with the conclusion of the kinetic 
formation model. However, the more realistic space-time evolution given by the transport model leads to an overall 
suppression of J /ip particles in the parton phase, instead of an enhancement predicted by the kinetic model. 

We have also shown that J/ip can be produced in the hadron phase at RHIC. This is due to the large D meson 
mass compared with charm quark mass. Although we have not taken into account parton-hadron interactions in the 
multi-phase model, both parton phase and hadron phase can exist in the same space-time region. This allows the 
production of J/ip^s in the mixed phase. At RHIC energies, the multi-phase model shows that hadron medium effect 
may lead to a net destruction of the J/ip produced in the parton phase. In addition, color screening is shown to have 
significant effect on the final J/ip yield. 

In the present study, we have not included the effect due to gluon shadowing p^ , which would reduce the charm 
quark and J/ip numbers. As the multi-phase transport model predicts an overall suppression of the directly produced 
J/ip, including gluon shadowing effect on the initial charm quark distribution will not change qualitatively the conclu- 
sion from the multi-phase model. Also, J/ip from decays of ip' or Xc was not considered in our study. The contribution 
from primary ip' or Xc to final J/ip are negligible due to their large destruction cross sections compared with that for 
J/ip. Although the production cross sections for ip' and Xc from charm quark interactions and charm meson interac- 
tions are expected to be larger than those for J/ip, the thresholds for these reactions are also higher. In particular, 
the statistical fragmentation model |^,^ indicates that the feeddown should be small for central collisions. We thus 
expect that our conclusion of J/ip suppression in heavy ion collisions at RHIC will not be qualitatively modified. 
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FIG. 1. Time evolution of the production rate (dashed line) and the destruction rate (dotted line) per unit rapidity for J/ip 
with \yj/^\ < 1 in the parton phase with the color screening effect. The thin dashed line denotes the collision rate per unit 
rapidity of charm and anti-charm quarks which would have produced the J/ip if color screening effect is absent. 
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FIG. 2. Time evolution of produced number (dashed line), destructed number (dotted line), and net number (solid line) per 
unit rapidity for J/tl) with \yj/^\ < 1 in the parton phase with the color screening effect. 




FIG. 3. Time evolution of the production rate (solid line) and the destruction rate (dashed line) per unit rapidity for J/ip 
with Ij/j/^I < 1 in the parton phase without color screening effect. 
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FIG. 4. Time evolution of produced number (dashed line), destructed number (dotted line), and net number (solid line) per 
unit rapidity for J/tp with \yj/^\ < 1 in the parton phase without color screening effect. 




FIG. 5. Time evolution of the production rate (solid line) and the destruction rate (dashed line) per unit rapidity for J/ip 
with \yj/^\ < 1 from the hadron phase for three values of D meson mass. Color screening effect is included in the initial parton 
phase. 
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FIG. 6. Time evolution of produced number (dashed line), destructed number (dotted line), and net number (solid line) per 
unit rapidity for J/ip with < 1 for three values of D meson mass. Color screening effect is included in the initial parton 

phase. 
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FIG. 7. Time evolution of the production rate (solid line) and the destruction rate (dashed line) per unit rapidity for J/'0 
with [j/j/^l < 1 from the hadron phase for three values of D meson mass. Color screening effect is not included in the initial 
parton phase. 
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FIG. 8. Time evolution of produced number (dashed line), destructed number (dotted line), and net number (solid line) per 
unit rapidity for J/tp with < 1 in the hadron phase for three values of D meson mass. Color screening effect is not 

included in the initial parton phase. 
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FIG. 9. Number of produced J/tp per unit rapidity with \yj/^\ < 1 as a function of D meson mass. Open symbols are for 
production in the parton phase while filled symbols are for the final number including production from the hadron phase. Circles 
and diamonds are, respectively, the results with and without color screening in the parton phase. Error bars are statistical. 
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FIG. 10. Midrapidity densities of charm quarks before (lower dashed hne) and after (lower solid line) parton evolution as 
well as of charm mesons before (upper dashed line) and after (upper solid line) hadron evolution. 
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FIG. 11. Midrapidity density of J /%!} as a function of the charm quark mass for the production in the parton phase (open 
diamonds) and the final number including the production in the hadron phase (filled diamonds). Error bars are statistical. 
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